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Abstract. I show that fractional exclusion statistics (FES) is manifested in general 
interacting systems and I calculate the exclusion statistics parameters. Most 
importantly, I show that the mutual exclusion statistics parameters-when the presence 
of particles in one Hilbert space influences the dimension of another Hilbert space-are 
proportional to the dimension of the Hilbert space on which they act. This result, 
although surprising and different from the usual way of understanding the FES, renders 
this statistics consistent and valid in the thermodynamic limit, in accordance with the 
conjucture introduced in J. Phys. A: Math. Theor. 40 F1013 (2007). 
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1. Introduction 

Fractional exclusion statistics (FES), introduced by Haldane in Ref. p] and with the 
thermodynamic properties calculated mainly by Isakov [2] and Wu [3], has received very 
much attention since its discovery and has been applied to many models of interacting 
systems (see for example Refs. |^ El El El El EIB IIIl 03 OS 0^ 03 OS IT3 OS 113 ED] ) - 
Several authors have also discussed the microscopic reason for the manifestation of FES 

CZIll22linilIElinilI5llIEllISlE3lCa3. 

Iguchi and Sutherland [23J showed that liquids of particles in three dimensions, 
interacting through long-range forces exchibit the nature of quantum liquids with FES, 
the characteristics of the FES being determined by the interaction. 

Murthy and Shankar [T7] analysed a system of fermions in the Colagero-Sutherland 
model. The system has a constant density of states (DOS) (along the single particle 
energy axis) and has a total energy of 

E({n l }) = J2e i n l + ^N(N-l), (1) 

where is the population of the single particle state of energy Cj, a = (e* — Cj_i) _1 
(for any i > 0) is the DOS, V is the mean-field interaction potential, and iV is the total 
number of particles in the system. By redistributing in an uneven way the interaction 
energy between the particles of the system and associating to the level i the quasi-particle 
energy 

e i = €i + Va- 1 J2n j , (2) 

3=0 

Murthy and Shankar obtained a gas with FES of parameter a — 1 + V. 

A model which is similar to that of Murthy and Shankar [17] has been employed 
also in Refs. [151 EE HH] to describe anyons on the lowest Landau level, coupled chiral 
particles on a circle, or interacting bosons in two-dimensions. 

In Refs. [24, 25J I showed that the same model, with a slight generalization, can 
lead to a condensation, which is a first order phase transition. 

In this paper I will extend the method of Murthy and Shankar to systems of general 
DOS and any interaction potential, (where i and j label the single particle states) and 
I will show that such systems lead to a more general manifestation of FES. While in the 
Murthy and Shankar model we have only direct exclusion statistics (i.e. the exclusion 
statistics is manifested only in the subspace where the particles are inserted) of constant 
parameter, a, here, in the general case, we shall have also mutual statistics (acting from 
one subspace into another); therefore we shall have more complex parameters, denoted 
as OLiy I will calculate explicitely the parameters o^- and I will prove that the mutual 
parameters (a^-, with i ^ j) are proportional to the dimension of the Hilbert subspace 
on which they act, verifying in this way the conjecture put forward in Ref. [26J. 
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Figure 1. The single particle energy levels in the noninteracting system (left) and the 
corresponding quasiparticle enegy levels (right, efj — + X)j=o ^i n i) when there are 
four particles on level 0, two on level 1, and one particle on each of the levels 2, 4, and 
5. In this particular case I chose 6j = i and Vij = 1/3 for any i, j. 



2. FES in systems with interaction 

Let us generalize the model of Murthy and Shankar [TTJ, [To] dSJ [19J by writing the total 
energy as 

E = J2 d n i + v ij n i n i (3) 

i ij 

and the quasiparticle energies as 

i-i x 

6 « = e i + V i3 n i + 7; V H n i ( 4 ) 

i=o ^ 

(see figured]). To make the calculations and the physical implications as clear as possible, 
we assume that we have bosons in the systems-in this way we shall not have to worry 
about adding a unit to the direct exclusion statistics parameters. I will also assume that 
the system is large enough, so that the spectrum is (quasi)continuous, of the (generally 
not constant) DOS, er(e). Then, assuming that Vy depends only on the energies of the 
interacting particles, in Eq. (j3j) I drop the subscript i and I use e as a variable, to write 

e = e+ /V(e, e>(e>(e') de 1 . (5) 

J 

In Eq. (j5J I also ingnored the term V(e, e)n(e). Although this term, for e = 0, may cause 
a first order phase transition [2U [25], here I just want to emphasize the characteristics 
of the emerging FES and earring along this term in the calculations would be useless. I 
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Figure 2. Left: e(e) before (solid curve) and after (dashed curve) the insertion of 
extra particles in the interval (ei,^)- This changes the values of ei, €2, and 63 into 
e'-i, e' 2 , and e 3 , respectively. To the right I draw the inverse of cr(e), to emphasize the 
change of the number of states (/ cr(e)de) in each of the intervals after the insertion 
of particles. Both plots are schematic and are used only to illustrate the principle of 
calculation. 

assume also that the function e(e) is bijective, so that I can use freely its inverse, e(e). 
Since e(e) and e(e) depend also on the populations of the energy levels below e or below 
e, respectively, I shall use also the notations e n ( e '< e )(e) and e n (e/<g)(e) whenever this will 
be needed for clarity. 

If I denote the density of states along the e axis by er(e) and the number of particles 
between the energy levels e\ and e%, by N(i\, £2), then we have the relation 

JV(ei,e 2 )= [ 2 a{e)n{e) de = f ^ a(e')n(e') de', 

Jei •'^(ei) 

where, obviously, n(e) = n[e(e)]. 

To show the underlying FES character of the system, I use the coarse- 
graining of the energy axis e. I split the quasiparticle energy axis into intervals- 
[eo, ei], . . . , [£m-i, ej\/], • • -which are small, but still contain large enough numbers of 
particles and energy levels; the FES will be manifested between and within these 
intervals [2H [22l HH HH HH HU HH H91 [231 120]. To each gj it corresponds an 
6j = e n ^ e<ei ){li) . I rewrite Eq. (jHJ) as a summation, 

M-l 

£m = e M + V(e M ,ei)N(ei,e i+1 ), (6) 

where, based on the fact that the intervals [ej,ej+i] are small and V is assumed to be 
continuous in both variables, I used the approximation V(eM, e») ~ V(€m, for any 
% < M. Using this decomposition I calculate the FES parameters. 

First I calculate the direct exclusion statistics parameter by adding Im-i particles 
in the interval [Im-i, £m]- Since I hold fix €m-i and €m, then €m-i and all the energy 
levels below it will also remain fix, while e« and all the energy levels above it will 
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change. I calculate the change of single particle states in the interval [em-x,cm\ ( see 
figure [21 with i = M — 1 and j = M); I use the notation e' M = eM,n(e'<e),i M -i(^)^ which 
is the value taken by em after the insertion of the Im-x particles. The initial number of 
states in the interval [e"m-x,cm] is G(em-x,^m) = J^f, ^K 6 ') ^' anc ^ after the addition 
of particles it changes into G'(em-x^m) — felLi a ( e ')de'. So, to calculate the difference 

SG(€m-x,€m) = G'(tM-x, cm) — G(em-x,£m) = Ism a { e ')& e ' , I calculate first the change 
of e' M : 

M-X 

e~M = 4/ + v ( e> M> £m-x)Im-x + J2 V ( e Mi £i)N(ei, e i+1 ) (7) 

i=0 

If I denote Sem = e' m — em and I expand V(em, around em, I get from (jHJ) and (JTj) an 
equation for 5e: 

,£ai-x)Im-x / x 
or, changing the summation into an integral, 

l+ d ^lM-x+f^ 9 -^*(e>)n(E')dE> [ J 

I look for linear effects, therefore I ignore the term proportional to Im-x from the 
denominator of equation (jUj) and I replace V(em, £m-i) by V(em, £m) (assuming that V 
is continuous in both variables). Writing 5G(em-x,^m) — • &(cm) = a e M eM^M-x, I 
get the direct exclusion statistics parameter 

om = Vi %' e ? )am ■ (10) 

Note that a u is identical to a calculated before [HI HSl HH H3 H51 HSJ if 9V[e, e(e,)]/de = 
0. 

Now let's calculate the mutual exclusion statistics parameters. For this I introduce 
Ii particles in the interval [ej,6j + i] (0 < % < M — 1). This will change all the energy 
levels Ej, of j > i (see figure [2]); let's denote the new values of Ej, j > i, by e'-. Taking 
all these into account, I write 

% M-X 

i = e' + V(e', U )U + £ VP, E 3 ]N(i 3l ~E j+1 ) + £ V[e', e>]N(~e„ e j+1 ). (11) 

3=0 j=i+X 

Expanding again V(e, e') to the linear order in both variables, I get the equation for 
= (4r)« - e M : 

,M-1 dV(e M ,ej) 



where I used the superscript to indicate that the particles were inserted at e*. The 
unknown quantities, 8e^ = (e'-)W — ej, can be calculated recursively, starting from 
j = i, using first equation (Q and then equation (1121) . By doing so, we first notice that 
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06j is proportional to Jj, for any j. Transforming both summations of equation ( 1121) 
into integrals and introducing the notation 

j{6M,£i) = j = j , l-Lo) 

I get the final equation for 5e, 

, ,~ ~ \ V(eM,ei) +/(eM,£i) r n ^ 

&(£M ' £i) = -TTjr~^w^ (14) 

If we plug in equation ffl3|) into equation (fT4l) . the later becomes an integral equation 
for <5e(e, e 

ng now the expression for St M , 
states in the interval [em-i, ^m\- 

do~(e 



Having now the expression for Sep}, we can calculate the change of the number of 



5G(e M -i, e M ) = cr(e M )5eM—cr(eM-i)$eM-i 



(^M-^M-i)Se M , (15) 



de 

CM 

where we ignored 8em — Sem, since 5 6m is itself a small quantity. Notice that because 
both, €m-i and em, vary at the insertion of particles at energies lower than em-i, the 
variation of the number of quasiparticle states in the interval [em-i, ?m] is proportional 
to Em — cm-1; i-e. is proportional to the dimension of the interval. Plugging equation 
( I14p into equation ffTBT) I obtain the mutual exclusion statistics parameter, 

dcr(e) 



(cm - e M -i){V(eM, gj) + fje, 
1 + Jo Hjr^ie'Me') 



(16) 



9e 

One can see immediately that if da(e)/de = for any e, as it was in the case of constant 
density spectrum, ai M i i = for any Im ^ h [HI [HI EEH 13 HSl [T9] . 

Now we observe directly the surprising character of the mutual exclusion statistics, 
namely that it is proportional to the energy interval on which it acts, (em — €m-i)- In 
Ref. [26] I showed that this characteristics is necessary to ensure the self-consistency 
of the FES formalism, especially in the thermodynamic limit. The method to calculate 
the particle population for such exclusion statistics parameters is also given there. 



3. Conclusions 



Fractional exclusion statistics (FES) is usually considered as an "exotic" type of 
statistics, manifested in special types of systems. Contrary to this belief in this paper, 
by analysing a system with a very general model of interaction between the constituent 
particles, I showed that FES is rather the rule than the exception. FES is manifested in 
general in interacting systems. Moreover, I calculated the FES parameters of the model 
gas and I showed that the mutual exclusion statistics parameters are proportional to 
the subspace on which they act. This conclusion is also in contradiction with the usual 
definition of FES and therefore seems peculiar. But it is not so. In Ref. [26] I showed 
that the typical definition of the mutual exclusion parameters leads to inconsistencies 
in the thermodynamics calculations and, in order to eliminate these inconsistencies, the 
exclusion parameters must have exactly the properties deduced here. 
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